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We study the Anderson localization in systems, in which transport channels with rather different 
properties are coupled together. This problem arises naturally in systems of hybrid particles, such as 
exciton-polaritons, where it is not obvious which transport channel dominates the coupled system. 
Here we address the question of whether the coupling between a strongly and a weakly disordered 
channel will result in localized (insulating) or delocalized (metallic) behavior. Complementing an 
earlier study in ID [H. Y. Xie, V. E. Kravtsov, and M. Miiller, Phys. Rev. B 86 014205 (2012)], 
the problem is solved here on a bilayer Bethe lattice with parametrically different parameters. The 
comparison with the analytical solution in ID shows that dimensionality plays a crucial role. In 
D — 1 localization is in general dominated by the dirtier channel, which sets the backscattering rate. 
In contrast, on the Bethe lattice a delocalized channel remains almost always delocalized, even when 
coupled to strongly localized channels. We conjecture that this phenomenology holds true for finite 
dimensions D > 2 as well. Possible implications for interacting many-body systems are discussed. 

PACS numbers: 72.15.Rn, 72.20.Ee, 72.70.+m, 73.20.Jc 



I. INTRODUCTION 

In a variety of physical contexts, the situation arises 
that two or more propagating channels with different 
transport properties are coupled together, competing 
with each other or modifying each other's properties. 
Under these circumstances it is interesting to study the 
resulting localization properties on the coupled system. 
Such a question arises in particular in the context of 
exciton-polaritons, which are hybrid particles: half pho- 
tons, half excitons, the two channels being coupled lin- 
early via dipolar interaction.^ Another realization of this 
physical situation can be found in bilayer graphene, with 
different degrees of disorder affecting the two layers. A 
recent work proposed such bi- or trilayers as field effect 
transistors, whereby a gate potential controls the degree 
of disorder sensed by the electrons in the bilayer.^ 

Similar questions arise in the problem of energy or 
matter localization in few- or even many-body problems, 
where a multitude of propagation channels may exist to 
transport particles or energy from one place in the sys- 
tem to another. For example, energy may be transported 
in small, nearly independent units in the form of quasi- 
particles, or it may have a propagation channel in which 
a larger amount of energy is propagating in the form of 
blobs of several quasiparticle-like excitations that form 
sorts of bound states. Such "bound states" were argued 
to be favorable transport channels in the context of few- 
particle problems. The problem was especially studied 
in low dimensions where under certain circumstances 
such compounds are found to have an enhanced local- 
ization length as compared to single-particle excitations. 
The question arises, then, as to which channel of propa- 
gation is the most favorable in transport problems con- 
taining a larger number of particles, or in the situation 
of particles at finite density. 

In this type of problem, the various propagation chan- 
nels are not independent, but couple to each other by 



scattering events. Understanding transport in such in- 
teracting systems is a challenging and largely unresolved 
problem. Here we do not aim at resolving all aspects of 
the many-body problem, but address one sub-question 
which arises in its context. Indeed, the interacting sys- 
tems have a common feature with noninteracting hybrid 
particles: Two or more propagating channels with para- 
metrically different localization properties are coupled to- 
gether and influence each other's transport characteris- 
tics. Under these circumstances it is interesting to study 
what are the resulting localization properties in the cou- 
pled system. In particular in the specific case where a 
less localized system is coupled to a more localized one, 
the question arises as to which of the two components 
eventually dominates the transport: Does one obtain an 
insulating or conducting system? A central result of our 
work is to show that the answer to this question depends 
crucially on the dimensionality of the system. 

In our recent worlP the question of the competition 
between alternative propagation channels was raised in 
disordered one-dimensional systems. This case can be 
studied in great detail in the form of a single-particle 
problem with two parallel, coupled channels. Among oth- 
ers, this naturally describes the Anderson localization of 
exciton-polaritons in quasi-one-dimensional semiconduc- 
tor heterostructures. By exactly solving the Anderson 
model on a two- leg ladder {D = 1)^ we found two regimes 
whose localization properties are qualitatively different: 
(i) a resonant regime, where the "slow" chain (the more 
disordered one) dominates the localization length of the 
ladder; this can be understood as a manifestation of the 
fact that in one dimension the backscattering rate de- 
termines the localization properties of a coupled system, 
since in general the localization length is of the order of 
the mean free path; (ii) an off-resonant regime, where the 
"faster" chain helps to delocalize the slow chain, although 
with low efficiency. 

In that ID study the disorder was taken to set the 
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FIG. 1. Schematic phase diagram for coupled Bethe lattices with identical hopping strength, but different disorder Wi / as 
inferred from the results in Figs. [S] and ^ The critical disorder for uncoupled lattices, Wc ~ 17.3, is indicated by the red lines, 
(a) Nearest-neighbor intralattice coupling only, 7 = 0. The mobility edge for the middle of the band {E = 0) of the coupled 
system is indicated by the black curve. In region A (yellow), in the absence of the coupling t± the two lattices would both 
be localized. However, the finite t± pushes the system into the delocalized phase. If in the absence of coupling one lattice is 
delocalized and the other one localized, there are two possibilities. In region B (gray), the coupled system becomes delocalized; 
that is, the less disordered channel dominates. Only when the delocalized lattice is very close to criticality and is coupled to a 
very strongly disordered lattice [region C (green)] localization prevails. However, this regime occurs in a very narrow window of 
parameters, (b) Next-nearest-neighbor hopping included, 7 = 1. The mobility edge is indicated by the blue curve. In region A 
(yellow) the coupling between two localized lattices induces a delocalized phase. In contrast to (a), the region C is eliminated 
by the next-nearest-neighbor hopping: The coupled system is always more delocalized than either of the uncoupled lattices. 



smallest energy scale, which allows for a fully analytic 
solution of the problem. In higher dimensions {D > 2), 
however, weak disorder has no significant effect on local- 
ization. Hence, we are restricted to considering interme- 
diate or strong disorder in order to address meaningfully 
the question of the role of interchannel coupling. Mean- 
while, since the disorder is comparable to or stronger than 
the hopping strengths, resonance conditions, as in regime 
(i) of the weakly disordered ID chains, are impossible. 
Furthermore, in contrast to the physics in one dimen- 
sion, proliferation of backscattering plays a subdominant 
role for the Anderson transition of the coupled system, 
and therefore, the resulting phenomenology of coupled- 



channel problems turns out to be rather different. 

In this paper we study two coupled Bethe lattices with 
different transport characteristics. This can be viewed as 
the limit of infinite dimensions {D 00) of the problem 
of coupled channels, which we will contrast with the case 
of two coupled chains {D = 1). The behavior on the 
Bethe lattice is suggestive of the physics to be expected 
in high-dimensional systems. Indeed, we believe that the 
qualitative behavior of coupled lattices in I) > 2 is very 
similar to the phenomenology found on the Bethe lattice. 
However, the latter has the significant advantage of being 
exactly solvable, which we exploit below. 



Statistical models on the Bethe lattic^^ have at- 
tracted a lot of studies, because they admit exact solu- 
tions and reflect features of the corresponding systems in 
sufficiently high spatial dimensions. The Anderson model 
on the Bethe lattice was first introduced and solved by 
Abou-Chacra, Anderson, and Thouless in Refs. [TOl and 
[TTJ where the existence of the localization transition was 



proven and the location of the mobility edge was found. 
That work showed in particular that localization is pos- 
sible in the absence of loops in the lattice. The model 
was solved by studying the self-consistency equation for 
the on-site self-energy, which leads to a nonlinear integral 
equation for the probability distribution function of that 
quantity. The transition from the localized phase to the 
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delocalized phase is characterized by the instabihty of the 
fixed point distribution of real self-energies with respect 
to a perturbation with infinitesimal imaginary parts of 
the self-energies. Physically, the latter describes an in- 
finitesimally weak coupling to a dissipative bath which 
allows for decay processes. The above instability signals 
that local excitations start coupling to a bath on sites 
infinitely far away, which signals their spatial delocaliza- 
tion. 

The stationary distribution function of the self-energy 
can be found numerically wit h the help of a population 
dynamics, or "pool," method.liSIlil The original work by 
Abou-Chacra et al. ^^l^has inspired a number of studies 
in both the physicJ^^HI^ and the mathematic^^^K^ com- 
munities. The recent work Ref. 16 points out that the 
Anderson model on the Bethe lattice may have a fur- 
ther transition within the delocalized phase and corre- 
sponds to a transition in the level statistics. Here, we 
focus however on the standard delocalization transition, 
as discussed by Abou-Chacra et al. 

In the present work we generalize the approach by 
Abou-Chacra et al. to the case of two coupled Bethe 
lattices. Following Refs. ITOl and [TT1 we derive a recursion 
relation for the local Green's functions (encoded in a 2 x 2 
matrix in layer space) and study the effect of interlayer 
coupling on the location of the transition. We restrict 
ourselves to the band center (£' = 0). Furthermore, we fo- 
cus on the case of lattices with identical hopping, but dif- 
ferent disorder strengths. This choice is motivated by the 
one-dimensional case, where equal hoppings lead to res- 
onance effects, which enhance the localization tendency 
in the coupled system. In contrast, we find that despite 
the choice of equal hoppings such a localizing effect al- 
most never occurs on coupled Bethe lattices. This is illus- 
trated by the schematic phase diagrams of Fig. [l] which 
anticipate and summarize the main results of our analy- 
sis: Under most circumstances the coupling between two 
layers enhances delocalization. Only when one couples 
a barely metallic layer to a strongly disordered second 
layer and excludes next-nearest-neighbor interlayer cou- 
plings [7 = in the Hamiltonian ([T]) below] , the coupling 
can induce localization. However, in the largest part of 
the phase diagram the coupling has a delocalizing effect. 
In the case of next-nearest-neighbor interlayer couplings 
(i.e., nearest-neighbor coupling across layers, 7 = 1), the 
coupled layers are always delocalized if one of the uncou- 
pled layers is delocalized. Moreover, in some range of 
parameters a coupling between two localized lattices can 
induce delocalization. 

Our central result may be summarized by the state- 
ment that on Bethe lattices the delocalized lattice essen- 
tially dominates the physics. In other words, if a delo- 
calized channel exists, delocalization, diffusion, and the 
ability of entropy production will persist even upon cou- 
pling to more localized channels. As mentioned before 
this is quite opposite to the phenomenology in ID where 
most often the more disordered chain dominates the lo- 
calization properties. 

The remainder of the paper is organized as follows. In 
Sec. [TT] we define the Anderson model on two coupled 




FIG. 2. Anderson model on a bilayer Bethe lattice, described 
by the Hamiltonian ([T]), shown for connectivity K-\-l — 3. We 
consider two types of interlayer coupling: (i) Only nearest- 
neighbor coupling (horizontal blue lines), setting 7 = 0. 
(ii) Additional next-nearest-neighbor coupling (diagonal green 
lines), with 7=1. 



Bethe lattices and derive the recursion relation for the 
local Green's functions. In Sec. HI we present the pop- 



ulation dynamics algorithm, which is used to study the 



statistics of the local self-energy. In Sec. IV we obtain the 
location of mobility edges, which gives rise to the phase 
diagrams shown in Fig. [l] Their qualitative features will 
be explained by a perturbative analysis. Finally, we dis- 
cuss the role of dimensionality and the possible impli- 
cations of our results on the interacting particles in the 
Conclusion. 



II. BILAYER ANDERSON MODEL 



Model 



We consider two Bethe lattices labeled by v = {1,2}. 
A Bethe lattice is defined as the interior of an infinite 
regular Cayley tree, each vertex having the same coor- 
dination number K -\- 1. The essential feature of such a 
lattice is the absence of loops. The Bethe lattice can be 
realized as the thermodynamic limit of a random regular 
graph of constant connectivity (i^T + 1), that is, a graph 
where each site connects to other K -\- 1 sites, which are 
randomly and uniformly selected. It is known that any 
finite portion of such a graph is a tree, with probability 
tending to one as the size tends to infinity. The advantage 
of the random-graph construction is the explicit absence 
of boundary effects. A random graph can thus be viewed 
as a regular tree wrapped onto itself. 

Analogously to the two-chain model studied in Ref. |3l 
we define the Anderson model on coupled Bethe lattices 
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as (cf. Fig. [2| 



(1) 



Here i labels the coordinates of two corresponding sites 
in the two layers, and (i, j) denotes two nearest neigh- 
bors i and j on the Bethe lattice. We take the onsite 
energies e^^^ to be independently distributed random vari- 
ables with zero mean, is the nearest-neighbor hopping 
strength within each layer. As motivated above, we take 
the intralayer hoppings to be equal, so as to come closest 
to the resonant case in one dimension, which shows the 
strongest localizing effects. t± is the interlayer hopping 
strength. In addition to direct (nearest-neighbor) inter- 
layer coupling, we also allow for next-nearest-neighbor 
hoppings of strength jt±. We will consider the two cases 
7 = and 7 = 1. For K = 1, the former reduces to the 
ID model studied in Ref. |3l 

The reason to include finite next-nearest-neighbor hop- 
pings is as follows. Consider the effect of coupling a first 
lattice to another one with very strong disorder or vanish- 
ing hopping. If we exclude next-nearest-neighbor hopping 
by setting 7 = 0, the only significant effect of the cou- 
pling is to increase the effective disorder on the first lat- 
tice, while the renormalization of its hopping is strongly 
subdominant or even absent altogether. However, the lat- 
ter is not the case if we allow for next-nearest-neighbor 
interlayer hopping. Indeed, this introduces a weak but 
nonnegligible additional propagation channel, which pro- 
ceeds via the disordered lattice. As the study below will 
show (and as anticipated in Fig.[T]), in the case 7 = 1, the 
renormalization of the hopping always dominates over 
the enhancement of effective disorder, and hence, the cou- 
pling always has a delocalizing tendency. In the context 
of the more general problems of coupled parallel propa- 
gation channels in many-body systems, the case of 7 7^ 
appears to be a rather generic and natural choice. Even 
a rather small 7 is sufficient to avoid the phenomenology 
found for 7 = 0, which leads to atypical behavior in some 
small regions of the phase diagram. 

In the Hamiltonian ([T]) the two layers are subject to 
different random potentials, characterized by two prob- 
ability distribution functions Pzy(e). For convenience we 
assume them to be box distributed: 



ee[-W,/2,W,/2], 
0, otherwise. 



(2) 



Our goal is to study the effect of weak interlayer cou- 
pling t± on the Anderson transition of the system. As 
mentioned above, this parallels the case of two resonant 
chains described in Fig. 7 of Ref. 3, However, as we will 
discuss in detail in Sec. |IV[ the parameter range of inter- 
est on the Bethe lattice is Wi^2 ^ t±it\\^ iii contrast to 
the weak disorder limit considered in Ref. |3l However, 



the notion of resonant interlayer coupling is meaningful 
only if the disorder is so weak that a well-defined disper- 
sion relation exists, which is not the case for the regime 
of interest on the Bethe lattice. Therefore, the equality 
of the two intralayer hoppings does not have important 
consequences in the present study. 



B. Recursion relation for the local Green's 
functions 



by 



The retarded Green's function at energy E is defined 



(3) 



where z^, /i G {1, 2} are layer labels, the kets stand for 

|z, i/) = cj^ I Vacuum), (4) 

and 77 is an infinitesimal positive real number, represent- 
ing an infinitesimally weak coupling to a dissipative bath 
into which particles can decay. We introduce 2x2 ma- 
trices in the layer space, Hi^ Gi^ and T, whose elements 
are 



and 



-^^ivt\\ - (1 - 5^v)lt^. 



(5a) 



(5b) 



(5c) 



T describes the hopping from one pair of sites to a neigh- 
bor pair of sites. 

One can easily show that Gi = Gi{E) satisfies the fol- 
lowing equation: 



G, 



1 



E^ir]-Hi-fYGff 



(6) 



jedi 



where di denotes the set of neighbors of i. The G^*^ are 
the Green's functions at the coordinate j in the absence 
of all bonds between the pairs of sites at i and j. G 
G^*^ {E) satisfies the recursion relation. 



E^IT]- 



(7) 



where dj \i denotes the set of neighbors of j excluding i. 
Gi and G^*^ are complex symmetric matrices. In order to 



obtain G^, we first solve the recursion relation ([7|), and 
then substitute the solution into Eq. (|6|. 

The self-energies are defined via the layer-diagonal ma- 



trix elements {G^^^)yy as 



Sj,{E) = E + ii^-e,,-l/{Gf) 



(8) 
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Their imaginary parts, 



(9) 



characterize the decay processes of local excitations over- 
lapping with I j, i^) and having energy E}^ 

Under the recursion ([7| the F^^ assume a stationary 
distribution, whose characteristics determine whether the 
system is in the localized phase or in the delocalized 
phase.^^ ^-*^*"*^^ * In the thermodynamic limit, one has 



lim lim P(ri > or > 0) 

77^0 AT^oo 



0, localized, 

> 0, delocalized, 



(10) 

JV being the number of lattice sites. The thermodynamic 
limit. A/" ^ 00, and the limit of vanishing dissipation, 
7^ ^ 0, do not commute, since in a finite system, whose 
spectrum is discrete, 77 ^ always leads to vanishing 
Tjjy^s. Note that the values r^^i^2 on the two sublattices 
are statistically dependent in the presence of coupling; in 
particular, they are of the same order of magnitude. 

We emphasize that the average value of Fj^, namely 
(Fj,), cannot be used to identify the Anderson transition, 
because in the localized phase an infinitesimal dissipation 
T] leads to long tails in the distribution function of F^^, 
which leads to a finite value (Fj^). Instead, one needs 
to consider the typical value of F^, as defined by the 
geometric average 



F, -e^^"^-> 



(11) 



which depends on the lattice label v if the two lattices 
are statistically not identical. However, as they are of 
the same order of magnitude, the localization transition 
can be identified by either of the two typical values, by 
locating the boundary between the two regimes: 



lim lim Ft,p. = r' (12) 
ri^^M^oo 1 > 0, delocalized. ^ ^ 



The equivalence of Ftyp,ry=i,2 for the purpose of iden- 
tifying the phase transition will be shown explicitly in 
Sec. |III[ based on the population dynamics. 



III. ANDERSON TRANSITION AND 
POPULATION DYNAMICS 

A convenient way to determine the mobility edge was 
proposed in Refs. [10] and [Til It is based on analyzing 
the stability of the real solution of Eq. ^ at the en- 
ergy E with respect to the insertion of the infinitesimal 
imaginary energy shift irj. In the localized phase the real 
solution is stable. In contrast, in the delocalized phase, 
the solution develops a finite imaginary part, which im- 
plies that Ftyp,ry(£^ irf) / as 7^ ^ 0. The physical 
interpretation of this criterion is as follows. For a fi- 
nite but large tree, if the boundary sites are coupled to a 
bath with a dissipation rate 7/, we test whether the dis- 
sipation at the root of the tree, measured by Ftyp,iy(£^), 



is vanishing or not as the tree size tends to infinity. If 
E belongs to the localized part of the spectrum (point 
spectrum), particle transport is absent at large scale and 
there is no dissipation at the root. In contrast, in the de- 
localized regime, we observe finite dissipation even deep 
inside the tree. This procedure in fact implements the 
criterion (10) for the Anderson transition, as the insta- 



bility of real self-energies refiects the Anderson transition 
as a phenomenon of spontaneous breakdown of unitar- 
ity of the scattering matrix associated with the system 
Hamiltonian.^^ 

The stability analysis can be realized by a popula- 
tion dynamics, which is a numerical recipe to solve the 
stochastic iteration Eq. A detailed description of 

such an algorithm for the single-lattice case can be found 
in Refs. 10 and 13HT61 The basic idea is to simulate 
the distribution of a random variable X by the empir- 
ical distribution of a large population of representatives 
{Xcy}-^^. Here the random variable X is the symmetric 



2x2 matrix G^^^ . For simplicity, we denote 0)^^ by G, and 

the population by {Gal^^. The M.^ 1 representatives 
can be understood as values of Green's functions on sites 
at a given distance from the root on a large tree. The 
population dynamics consists of a number of sweeps of 
the population, which simulate the propagation of dissi- 
pation step by step towards the root of the tree, whereby 
the number of representatives is kept constant. At the 
71*^ stage, we denote the population as {60,^^^}^! 7 which 
are obtained with the following procedure: 

(i) As an initial condition for the population we chose 
the Green's functions of M uncoupled sites subject to 
a random potential and a small dissipation, that is, 
{Gq,,o}^i with matrix elements 

(Gc,,o).. = {E- e,, + ir])-\ v e {1, 2}, (13) 

(Ga,o)l2 = (Gq;,o)21 = 0, 

where e^u are independently drawn from the probability 
distribution ([2|. 7^ is taken as small positive number, 
representing the dissipation on the boundary sites of the 
tree. 

(ii) Generate the nf^ population from the {ug — 1)*^ 
population. For each member P = 1,2,-- - of 
the new population, we choose K matrices randomly 
and uniformly from the population {Ga,ns-i}'^i^ called 
{Gci,n,-i, • • • ,Gc^,n,-i}, and generate 2K random 
numbers according to the probability distribution func- 
tion in Eq. ([2]), called {ei^^, • • • , cku} with u = 1,2. Sub- 
stitute these quantities on the right-hand side of Eq. ^ 
with 7^ = since the dissipative bath only couples to the 
boundary sites, and obtain G/3,ns on the left-hand side. 

We calculate the typical value of Tj^ in the population 

{Ga,ns}c 



^(0 



InF 



1 

M 



M 



ElnF^"" 



(14) 



using Eqs. ([8]) and ([9]). The localization transition can 
be determined by studying the evolution of F^yp^^ un- 
der sweeps. As shown in Ref. H^O^, if the Vj^s are small 
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enough, the recursion relation ([7| leads to a linear homo- 
geneous equation for Fjjy, and the growth of the typical 
value of Tjiy under sweeps is dominated by the largest 
eigenvalue of the linearized recursion relation. Therefore, 
as long as ^typ]y is sufficiently small, statistically ^typ]u 
grows linearly with the growth rate 



InT 



ins) 



inr 



{ris-l 



(15) 



Notice that in this linear regime as long as the two lattices 
are coupled, the statistics of A^^ is independent of the lat- 

(n ) 

tice index u. In other words, T^yp^^^^^ 2 deviate from zero 
simultaneously as the system crosses into the delocalized 
phase, and therefore the criterion for delocalization tran- 
sition ( p!2| ) does not depend on u. The statistical analysis 
of Xn^ below is restricted to the linear regime where the 
remain small. 

The average growth rate of F^^p'^^ over ?ig ^ 1 succes- 
sive sweeps is given by 



(16) 



and the standard deviation is 



5\ 



(17) 
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Physically, |A|~^ may be interpreted as a localization 
length in the insulating phase, or as a correlation length 
in the delocalized phase. The Anderson transition occurs 
wheiP^ 



A = 0. 



(18) 



We obtained numerical results using a population size 
M = 10^, dissipation r] = 10~^^, and = 200 sweeps. 
The statistics of A^^ was collected only after about 10 
sweeps to avoid the initial transient. We checked that 
the 77 dependence of A and 5X was very weak, as long as 
?7 was taken to be small enough. 



FIG. 3. Numerical results for the growth rates A [Eq. (16)] 
at the band center for statistically identical Bethe lattices as 
functions of disorder strength Wi = W2 ■ Energies are in units 
of t|| = 1. The error bars correspond to ^A [Eq. (17)]. For the 
uncoupled lattices (red triangles) the critical disorder strength 
is Wc{t± = 0) ^ 17.3. Upon coupling the lattices we find the 
critical disorder strengths: (a) Wc{t± = 1,7 = 0) ~ 20.7 
(black squares); (b) Wc{t± = 1,7 = 1) 37.5 (blue circles). 



A. Numerical results 

In the numerical calculations we took a connectivity 
K -\- 1 = 3 and hopping strengths t± = = 1. We 
analyze the two cases in turn. 



IV. PHASE DIAGRAM 

Let us now analyze the effect of the interlayer coupling 
t± on the Anderson transition. For convenience we focus 
on the band center, E = 0. We concentrate on rela- 
tively weak interlayer coupling tj_ < ty, which guarantees 
that in the absence of disorder the energy bands are not 
substantially changed by the coupling. In this case a mo- 
bility edge first appears at the band center E = upon 
increasing the hopping strength.!^ Below we present the 
numerical results of the population dynamics, which lead 
to the phase diagrams shown in Fig. [l] In Sec. |IVB| a 
perturbative analysis is given to explain the qualitative 
features of the phase diagram. 



1. Statistically identical lattices (Wi = W2) 

We first analyze two statistically identical lattices with 
disorder strength W = Wi = W2 (following the di- 
agonal line in Fig. [I]). In Fig. [3] we show X ± SX at 
£^ = as functions of the disorder strength for uncou- 
pled and coupled lattices. The transition point is de- 
termined by Eq. (18). For the uncoupled lattices we 

17.3, which 
For cou- 



find the critical disorder Wc{t± = 0) 
agrees with the results in Refs. 15 and 16, 
pled lattices, the critical disorder strength increases to 
Wc{t± = 1,7 = 0) ~ 20.7 with nearest-neighbor coupling 
only, and to Wc{t± = 1, 7 = 1) ^ 37.5 when next-nearest- 
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neighbor coupling is included. Thus the critical disorder 
is enhanced by the coupling, Wc{t± 7^ 0) > Wc{t± = 0). 
This implies that if two decoupled lattices are in the lo- 
calized phase but close enough to criticality, the coupling 
will delocalize the system. 



2. Parametrically different lattices (Wi / W2 ) 

Let us now study two Bethe lattices with identical hop- 
ping but different disorder strengths. We take Wi = 
Wc{t±_ = 0) ~ 17.3 to be critical (following the red line 
Wi = Wc in Fig. [T]) and analyze whether the coupling 
to a more disordered lattice pushes the system to a lo- 
calized or delocalized phase. If the inter layer coupling is 
weak < for both 7 = and 7 = 1, we expect that 
the system is delocalized when W2 is not much larger 
than Wc- This is expected from the results of the pre- 
ceding section. However, for W2 ^ Wc the situation may 
depend on the type of interlayer coupling. 

In Fig. [ij we show A =b as functions of W2- We 
observe the following features: For 7 = a mobility edge 
occurs at the fairly large disorder W2 = T^2,c(^± — ^-,1 = 
0) ^ 47. In other words, as long as W2 < W2^c the band 
center becomes delocalized, while it is localized beyond 
VF2,c- However, when next-nearest-neighbor hopping is 
included, with relative strength 7 = 1, the band center 
becomes always delocalized upon coupling, for any value 
of W2' As W2 00, the two lattices decouple effectively, 
and the band center tends back to criticality, from the 
localized and the delocalized side, for 7 = and 7 = 
1, respectively. Empirically we find that A ~ c{j)/W2 
for large W2, where 0(7 = 0) < and 0(7 = 1) > 0. 
As will become clear from the perturbative analysis in 
Sec. |IVB[ this is due to the suppression or enhancement 
of the probability of resonances between two neighboring 
sites on the first lattice. That effect is of the order of 
I/W2. 

The results obtained in Figs. [3] and H| give rise to the 
schematic phase diagram shown in Fig.|lJ One can distin- 
guish three regions according to the effect of the interlayer 
coupling: 

(i) Region A (yellow area). In the absence of coupling 
the two lattices are both localized but close enough to 
criticality. The coupling pushes the two nearly critical 
lattices into the delocalized phase. 

(ii) Region B (gray area). The better conducting lat- 
tice is (sufficiently far) in the delocalized phase, while the 
more disordered lattice is strongly localized. The coupled 
system is nevertheless delocalized due to the dominance 
of the better channel. 

(iii) Region C (green area). In the absence of coupling 
the less disordered lattice is delocalized but very close to 
criticality. The more disordered lattice is strongly local- 
ized. If there is nearest-neighbor coupling only (7 = 0), 
it pushes the system to the localized phase. However, 
this atypical region is entirely absent if a strong enough 
next-nearest-neighbor coupling is included (7 = 1). 



1^ 
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0.05 



0.00 
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1.95 2. 2.05 
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W2c~47.0 
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(b) 7 = 1 tii=t^=l Wi = W, 




FIG. 4. Statistically nonidentical Bethe lattices: Numeri- 
cal results for the growth rate A [Eq. (16)] at the band cen- 
ter as functions of the disorder strength W2- The disorder 
strength on the 1 lattice is fixed at Wi = Wc{t± = 0) ^ 17.3 
and the interlayer coupling is = 1. The other parame- 
ters are the same as in Fig. [3] (a) For the nearest-neighbor 
coupling there is a mobility edge, W2,c(t± = 1,7 = 0) ~ 47. 
As W2 00 the system approaches criticality from the lo- 
calized phase, and A ^ —I/W2, as expected analytically, (b) 
With next-nearest-neighbor coupling the system is always in 
the delocalized phase and approaches the transition point like 
A ^ I/W2 as W2 00 (best fit shown as dashed line in the 
log-log plot). 



B. Perturbative analysis 

The salient features of the phase diagrams shown in 
Fig. [1] can be understood qualitatively by applying a per- 
turbative analysis in the limit W2 ^ t\\,t±. The cou- 
pling to the strongly disordered second lattice has two 
competing effects on the first lattice: On the one hand, 
the hopping strength t\\ is effectively enhanced. On the 
other hand, the variance of the on-site energies on the 
first lattice is effectively enhanced, too. If the relative en- 
hancement of the hopping dominates, the coupling tends 
to delocalize the system. 

To leading order in I/W2, the correction for the hop- 
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ping strength between nearest-neighbor sites \i, 1) and 
|j, 1) is 



W2 



(19) 



where 



E-ti2 E-€ 



i2 



(20) 



Likewise, the correction of the local potential on site |i, 1) 
due to self-energy effects is 



t 



(21) 



where 



W2 W2 
^'''^^^2' = ^^^^ 

^1,2,3,4 are dimensionless random varia bles whose proba- 
bility distributions have long tailsP^ Hence both ^t\^i^ 
and btix are dominated by rare, large values. This im- 
plies that the dominant events are those where either 
the hopping strength or the disorder strength is strongly 
enhanced, that is, a link is either strongly favored or 
blocked. The ratio of the probabilities of such enhance- 
ments determines which effect is dominant. Notice that 
a typical value of ^tx^ij is of order 0(1/1^1 ) for 7 = 0, 
but 0(1/1^2) for 7 > 0, while Je^i scales as 0(1/ W2) re- 
gardless of the value of 7. Therefore, the enhancement 
of disorder is dominant when 7 = 0, that is, when the 
nearest-neighbor interchain hopping is suppressed. 

Let us now discuss the Anderson transition at £^ = 0. 
We base this discussion on two observations: First, the 
delocalization of wavefunctions on the Bethe lattice has 
recently been shown to occur along sing le paths.^ We 
should therefore study the decay rate of excitations along 
the best possible path for propagation. To obtain a qual- 
itative understanding of the effects of coupling, we ap- 
proximate the propagation amplitude between two re- 
mote sites of the first lattice as the product 



n 



^11 +^^i,ij(7) 



L > 1, (23) 



where L is the distance between the two sites. ni=i ^|| / ^ii 
is the amphtude in the absence of couphng, and 



^ ^ -rj l + '^^i,»j(7)A|| 
11 l + 5eii(7)/eii 



(24) 



represents the enhancement due to the coupling to the 
second lattice. The amplitude is the lowest order term 
in an expansion in the hopping. It corresponds to Ander- 
son's "upper limit" approximation, which neglects self- 
energy effects from sites lateral to the considered paths. 



as well as the regularization of resonances due to higher 
order corrections from hoppings along the path. Based on 
this approximation one obtains a simple approximate cri- 
terion for localization: Consider the probability P\Ai^\>\ 
that the propagation amplitude \Ai\ along the most fa- 
vorable path exceeds some fixed finite value 0(1). Local- 
ization obtains so long as this probability vanishes in the 
thermodynamic limit, P\Ai^\>i ^ as L ^ 00.^ 

In order to understand the phase diagram in Fig. [ij let 
us consider critical disorder on the first lattice, W\ = Wc^ 
and study how the probability P\Ae\>\^2) depends on 
W2 via the correction factor Rj^. For W2 ^ 00, Rl ^ I 
(in probability), and P\Al\>i(^'^ ~^ ^) behaves criti- 
cally; that is, it does not decay exponentially with L. 
For finite but large 1^2, we need to estimate the correc- 
tion factor Rl. As mentioned above, it is dominated by 
the rare events in which either the hopping strength or 
the local disorder are strongly enhanced, such that one 
of the factors in Eq. (24) is significantly different from 
unity. 

The probability of a strong enhancement of hopping, 
^\stl,^J\>tV ^^^^^^ ^^^^ ^2 like^"^ 



\Sti,ij\>t\\ 



\nW2 
Wf 



7 = 



7^ I 1 



(25) 



where co,i = 0(1) do not depend crucially on the pa- 
rameters of the system. The probability of the strong 
enhancement of local disorder behaves like 



^2(7) 



(26) 



where 02(7) = 0(1) depends on 7, and has a finite limit 
C2(7^0) >0. 

Among the L factors of Rl a fraction of order 
^dti ij\>t\\ is significantly larger than unity, and a frac- 
tion of order P\deii\>\eii\ terms significantly smaller than 
unity. Therefore, it is reasonable to assume that a typical 
value of Rl takes the form 



where the Lyapunov exponent 5(7, W2) is 



(27) 



(28) 



with of order 0(1). Substituting Eqs. (25) and (26) 
in Eq. (28), we predict the scaling 



sh,W2) 

with a coefficient 

C(7) r 



W2 ' 



W2 > Wi,tn 



/(7) 



Wi 



(29) 



(30) 



where f{j) (x 7 for 7 <C 1. Obviously, the condition 
0(7) = marks the transition between enhanced and sup- 
pressed propagation. Close to that criticality, the inverse 
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localization or correlation length follows from the growth 
rate |A| [cf. Eq. (16)] which is proportional to s(7, W2), 



A cx 5(7, W2) 



W2 



(31) 



The scahng with 1/1^2 is clearly observed in the numer- 
ical data of Fig. [4j confirming the dominance of rare 
events. 

Let us now discuss the 7 dependence of 0(7). Without 
next-nearest-neighbor interlayer hopping, we have c{j = 
0) ^ —7.5 < 0, as we numerically obtain in Fig. Qa). 
This is due to the fact that P^^t^ ij\>t\\ parametrically 
smaller than P\5eii\>\eii\ ^2 ^ 00 [cf. Eqs. (|25| and 
(26)]. Therefore, for large enough 1^2 (> ^2,c) the more 



disordered lattice drives a less disordered, critical lattice 
to the localized phase, as seen in regime C of Fig. [l] (a). 

However, when 7 > 0, the probabilities for signifi- 
cant corrections Sti^ij and ^e^i both scale as I/W2. For 
large enough 7, 0(7) becomes positive, as one may an- 
ticipate from Eq. ([30|, considering that ty/Wi is nu- 
merically small at criticality. Indeed, the case 7 = 1 
shown in Fig. [ijb) is already deep in this regime, with 
c(7 = 1) ^ 12.8 > [cf. Fig.gb)]. 

The equation 0(7 = 7c) =0 has a solution for some 
< 7c < 1. 7c determines the minimal next-nearest- 
neighbor interlayer hopping which assures delocalization 
upon coupling to a disordered lattice even in the limit 
W2 ^ Wi. A naive linear interpolation between 0(7 = 0) 
and c(7 = 1) allows us to obtain a rough estimate 



7c ^ 0.37 



(32) 



for the Bethe lattices of connectivity K-\-l = 3 considered 
here. 



V. DISCUSSION AND CONCLUSION 

We have studied the Anderson localization problem 
on two coupled Bethe lattices, which represents a two- 
channel problem in the limit of infinite dimensions. Our 
main result is the finding that a conducting transport 
channel is hardly ever localized by the coupling to more 
disordered channels. Rather, transport is usually en- 
hanced by such a coupling. This holds true except in the 
case where three conditions are met simultaneously: (i) 
the conducting channel is very close to criticality; (ii) it is 
coupled to a strongly localized channel; (iii) next-nearest- 
neighbor interlayer couplings are strongly suppressed or 
absent. Only in these exceptional cases the coupling to 
localized channels may induce a localized phase in an oth- 
erwise conducting channel. The coupling between mod- 
erately localized channels may instead induce delocaliza- 
tion. We believe that these trends persist also in high 
but finite dimensions {D > 2) where the metal-insulator 
transition takes places at strong disorder. This conjec- 
ture is based on the observation that in higher dimen- 
sions, as well as on the Bethe lattice, delocalization is 
mostly driven by a sufficiently strong forward scattering. 



whereas weak localization effects and enhanced backscat- 
tering play a much less important role than in D < 2. 
We believe that this difference is at the root of the very 
different phenomenology between coupled Bethe lattices 
and ID chains. 

In two dimensions, the localization length becomes 
parametrically larger than the mean-free path at weak 
disorder. However, since the proliferation of weak local- 
ization and backscattering leads to complete localization 
(in the absence of special symmetries), we expect that 
a well propagating channel becomes more strongly local- 
ized upon resonant coupling to a more disordered chan- 
nel, similarly as in one dimension. It might be interesting 
to investigate this numerically. Apart from its theoretical 
interest, the physics of coupled, unequally disordered 2D 
lattices might also have practical applications. For ex- 
ample, it was recently proposecP that a sheet of bilayer 
graphene with different disorder strength on the two lay- 
ers could be operated as a field effect transistor, whereby 
a perpendicular gate bias tunes the effective disorder of 
carriers. 

The study of localization properties of few- or many- 
particle systems is more subtle than the toy problem 
which it motivated here in part. The reason is that mul- 
tiple (much more than two) coupled channels with com- 
plicated substructures may exist to transport particles or 
energy. In particular, for few-particle problems, it has 
been shown that some of the channels, in which a large 
number of quasiparticle-like excitations are propagating 
in the form of "bound states," can be more efficient for 
transport than others, in which the excitations propa- 
gate essentially as independent units.^^ This suggests 
that one might have to think of the many-body prob- 
lem as having a hierarchy of channels with parametri- 
cally different transport properties. It is reasonable to 
assume that the effective dimensionality of such channels 
should be the same as that of the system. Our analysis 
of a two-channel model in ID has demonstrated that a 
"bad" channel dominates only when it is resonantly cou- 
pled to a better ("faster") channel. If the two channels 
are far from resonance, or if they live in higher effective 
dimensionality, the fast channel almost always dominates 
the localization properties, as the present study suggests. 
Moreover, even in ID resonance conditions are not met 
very easily. It either requires two channels with equal 
hopping strength, or an energy close to the band center 
or the band edges. Summarizing these considerations, we 
come to the qualitative conclusion that, apart from some 
exceptional cases, better conducting channels generically 
dominate the delocalization: A diffusing channel is diffi- 
cult to shut down by coupling to dirtier channels. 

In the context of interacting many-particle systems the 
above leads to the following conjecture: In order to estab- 
lish that a many-particle system conducts and is not fully 
localized, a sufficient condition will be found by identi- 
fying the best transport channel and showing that it is 
delocalized. Indeed, our study suggests that the inclu- 
sion of coupling to other channels usually only enhances 
transport. This observation should be a central ingredi- 
ent when generalizing the ideas of Refs. |4] and [5] to the 
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analysis of quantum dynamics and transport of systems 
with several particles. However, at this stage, the ap- 
plication to many-particle systems remains a conjecture 
which needs to be tested further. For example, one should 
establish whether coupling to a much larger number of 
slow channels does not alter our qualitative findings of 



"the survival of the fastest" . 
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